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AERODYNAMIC CHARACTERISTICS AND FLYING QUALITIES OF A 

TAILLESS TRIANGULAR— WING AIRPLANE CONFIGURATION AS 

OBTAINED FROM FLIGHTS OF ROCKET-PROPELLED MODELS 

AT TRANSONIC AND LOW SUPERSONIC SPEEDS 

By Grady L. Mitcham, Joseph E. Stevens, 
and Harry P. Norris 

SUMMARY 


A flight investigation has teen made at the Langley Pilotless 
Aircraft Research Station, Wallops Island, Va., to determine the aero- 
dynamic characteristics of models of a -tailless, triangular^wing air- 
plane configuration. The results from three successful flight tests are 
presented • for the Mach number range Between' 0.75 and 1.28. 

The data showed that the models tendedL to tuck under slightly 
through the transonic region. The variation of lift coefficient with 
angle of attack was linear within the range of angles tested and the 
lift— curve slope increased gradually between Mach numbers 0.88 and 1.00. 

The hinge-moment coefficients increased severely between Mach num- 
bers O.85 and 1.15 but showed a gradual decrease above a Mach number 
of 1.20. Elevator effectiveness decreased approximately 40 percent 
through the transonic region. 

The models exhibited static and dynamic longitudinal stability 
throughout the test Mach number range with the center of gravity located 
at 20 and 25 percent mean aerodynamic chord. The aerodynamic center 
showed a gradual rearward movement of about 15 percent mean aerodynamic 
chord in the transonic region. 

All the models possessed directional stability throughout the 
angle— of-attack and speed ranges of the flight tests. 
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An analysis of the flying qualities of a full-scale configuration 
has been made from the data obtained from the three flight— test models. 
The analysis indicates adequate elevator control for trim in level 
flight over the speed range investigated. Through the transonic range 
there is a mild trim change with a slight tuck Ing-under tendency. The 
elevat on-control effectiveness in the supersonic range is reduced to 
about one-half the subsonic value although sufficient control for maneu- 
vering is available as indicated by the fact that 10° elevator deflec- 
tion would produce 5g normal acceleration at a Mach number of 1.2 
at h0, 000- foot altitude. The elevator control forces are high and indi- 
cate the need of a control— boost system as well as the power required 
of such a system. The damping of the short— period oscillation is 
adequate at sea level but is reduced at 4-0,000 feet. 


INTE0DUCTI0N 


The NACA is conducting a flight investigation of rocket— powered 
models of a triangulai^jlng, tailless airplane configuration to 
evaluate stability and control at low supersonic and transonic speeds. 
Eesults obtained from the successful flight tests of three models 
served as a basis for the analyses presented in this paper. 

The models had a wing of triangular plan form with 60° sweepback 
of the leading edge and an aspect ratio of 2.31, the profile at al 1 
spanwise stations being an NACA 65 (q£)- 006.5 section. Longitudinal 

control was provided by a single set of constant-chord trailing-edge 
control surfaces on the wing called elevons. Deflecting the elevons 
together provides longitudinal control and, in the hypothetical air- 
plane, deflecting them differentially would give lateral control. The 
vertical fin of the models was of triangular plan form with a leading- 
edge sweepback of 60°. 

For the flight tests of the models, the program for control move- 
ment included abrupt pull-ups and push— downs with the elevons operated 
as elevators. The present paper contains the results of an analysis 
of the aerodynamic characteristics and the stability derivatives 
evaluated from the flight tests of the three rocket-powered models and 
an analysis of the flying qualities of such an airplane in the Mach 
number range from 0.75 to 1.2. The flying qualities are based on a 
hypothetical triangular-^wing airplane with a wing loading of 27.3 pounds 
per square foot at sea level and 40, 000-foot altitude. The computations 
are based on two center-of-gravity positions, 20 and 25 percent of the 
mean aerodynamic chord. 
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SYMBOLS 
time from launching, seconda 
Eeynolds number 

velocity, feet per second 
velocity of sound, feet per second 
Mach number 

free-etream static pressure, pounds per square foot 

specific heat ratio (value taken as 1.4o) 

hinge moment, inch— pounds; total impact pressure 
(equations (l) and (2)), pounds per square foot 

stick force, pounds 

stick movement, inches 

dyn am ic pressure, pounds per square foot 

mass density of air, slugs per cubic foot 
weight of model, pounds 
wing area (6.25 sq ft) 
mean aerodynamic chord (2.19 ft) 

longitudinal acceleration, feet per second per second 
normal acceleration, feet per second per second 
transverse acceleration, feet per second per seco nd 

acceleration due to gravity (32.2 ft/sec^) 
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angle of attack as measured during flight, degrees 

angle of attack corrected to center-of— gravity position, 
degrees 

distance "between center of pressure of ■ angle— of-attack 
vane and center of gravity of model, feet 

control deflection measured on chord line parallel to 
the plane of symmetry (positive with trailing edge • 
down), degrees 

chord— force coefficient, positive in a forward 
direction 


(H W l\ 
S q_) 


normal— force coefficient 


/fn W l\ 

\g S qj 


lift coefficient (°N cos a + C c sin 


pitching-moment coefficient 


a) 


rate of change of lift coefficient with angle 'of attack, 
per degree 

rate of change of lift coefficient with elevator deflec- 
tion for a constant angle of attack^ per degree 

trim angle of attack, degrees 
trim elevator deflection, degrees 


trim lift coefficient 


rate of change of total lift coefficient "between two trim 
conditions or elevator deflections, per degree 

rate of change of hinge-moment coefficient with angle 
of attack, per degree 

rate of change of hinge-moment coefficient with elevator 
deflection, per degree 
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■basic hinge-moment coefficient at zero angle of attack 
and zero elevator deflection 

rate of change of angle of attack with elevator deflect icn 
between two trim conditions 

rate of change of pitching-moment coefficient with angle 
of attack, per degree 

rate of change of pitching-moment coefficient with elevator 
deflection for constant angle of attack, per degree 

basic untrimmed pitching-moment coefficient at zero angle 
of attack and zero elevator deflection 

moment of inertia about pitch axis, slug-feet^ 

period of an oscillation, seconds 

time to damp to one-half amplitude, seconds 

cycles for the short-period oscillation to damp to 
one— tenth amplitude 

amplitude of a short— period oscillation 
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rate of change of pitching-moment coefficient with lift 
coefficient 


rate of change of elevator deflection with angle of attack 
(due to flexibility of control system) 

mass of model, pound— seconds^ per foot 
, per radian 


as used in equation ( 5 ), seconds per degree 
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ddc/2V 


per radian 
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as used in equation (5), seconds per degree 
pitch angle, degrees 

nondimensional angular velocity of pitch 
nondimensional rate of change of angle of attack 

Subscripts : 

a full-scale airplane 

m . model 

The symbols . and . . over a quantity represent, respectively, 
the first and second derivatives of the quantity with respect to time. 
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MODELS AND APPABATUS 
MODELS 


A three— view drawing of the models used in the present investi- 
gation is given in figure 1 and the physical characteristics of the 
models and a full-scale representative tailless triangular-wing air- 
plane are presented in table I. Photographs of one of the models are 
shown in figures 2 and 3. The model fuselage and components were con- 
structed of duralumin, magnesium castings, and magnesium skin. The 
fuselage construction was of the monocoque type divided into three 
sections. The three sections were the nose section which held the 
telemeters, the center section which held the wings, vertical fin, 
compressed-air supply, and control-actuating system, and the tail sec- 
tion which contained the rocket motor and booster attachment. 

The planned movement of the elevons called for abrupt pull-ups and 
push-downs operating at a frequency of about one cycle in 1.2 seconds 
and was accomplished by a compressed-air system. The control surfaces, 
which were unsealed, moved together between stops in an approx imat ely 
square-wave motion. On model 1 the surfaces were deflected down 5.3° 
•and up 5. 3°; on model 2 the deflection was down 4.7° and up 4.7°j and 
on model 3 the deflection was down 1.1° and up 5*2°. The controls were 
in operation during the entire flight. 
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The models were "boosted to supersonic speeds by a solid— fuel, 

6— inch-diameter Deacon rocket motor which is capable of producing an 
average thrust of 65OO pounds for approximately 3*1 seconds. 

The rocket— sustainer motor for the model was a 5-inch solid- 
fuel HVAB shortened to 17 inches and modified to give an average thrust 
of 900 pounds for 1.4 seconds. The small sustainer motor served a two- 
fold purpose; during the power-on portion of the flight it prevented 
immediate deceleration after . separation, allowing the controls to operate 
one complete cycle at approximately a constant Mach number, and it assured 
a positive separation between model and booster at booster burnout. 

The sustainer motor nozzle served as the point of attachment of the 
booster to the model. This type of attachment also allowed a separation 
of the booster from the model if the ratio of drag to weight of the model 
and booster were favorable. 

The booster-model combination was ground launched from a crutch- 
type launcher, as shown in figure 4. The launching angle from the hori- 
zontal for model 1 was 43°40* , for model 2 was 44°40' , and for model 3 
was 43°23* . Table II presents the weight and balance data for the models 
and the full-scale airplane. Figure 5 shows a sequence of photographs 
of one booster-model combination at take-off. 


APPARATUS 


The data from the flights were obtained by the use of telemeters, 
Doppler velocimeter radar, photography, and radiosondes. The time his- 
tories of the data as the models traversed the Mach number range were 
transmitted and recorded by a telemeter system which gave eight channels 
of information. The data recorded were longitudinal, transverse, and 
normal acceleration; hinge moment; control position; angle of attack; 
total pressure, and a reference static pressure used to determine free— 
stream static pressure. Figure 6 shows the instrumentation arrangement 
on a typical model. Angles of attack were obtained by a vane— type angle— 
of -attack indicator located on a sting ahead of the nose of the model. 

A description of this indicator can be found in reference 1. The angle— 
of -attack range covered by the indicator with the vane located on the 
center line of the model was approximately ±15°. On model 3 "the angle— 
of -attack sting was deflected down 10° from the center line of the model 
in order to record higher positive values of angle of attack. Figure 7 
shows a photograph of model 3 equipped with the deflected sting. 

Fixed wide-angle cameras and l6-millimeter motion— picture cameras 
recorded the launchings. The motion— picture cameras also tracked the 
flights. 
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TEST TECHNIQUE 


The models were disturbed in pitch by the abrupt movement of 
elevons operated as elevators at preset time intervals which gave an • 
approximately square-wave type of elevator motion. The desired aero- 
dynamic coefficients and longitudinal— stability derivatives were obtained 
by analysis of the hinge moments, angle of attack, and acceleration 
responses resulting from these cyclic disturbances. 


BASIS OF ANALYSIS 


The aerodynamic coefficients, stability derivatives, and flying 
qualities presented in this paper were reduced. from the model— flight 
data. The most recent specifications for satisfactory flying qualities 
(references 2 and 3) have been used in this analysis of flying qualities. 
Inasmuch as these specifications are restricted to subsonic speeds and 
the current range of interest is in the transonic speed range, no 
detailed step— by— step comparison with these specifications has been 
attempted. Figures 8, 9, and 10 are typical portions of the recorded 
time histories at low supersonic, transonic, and high subsonic veloci- 
ties. The variations in stability and control effectiveness through 
the Mach number test range can be seen in these figures by comparing 
the periods, and amplitudes of the short— period longitudinal oscillations. 
A discussion of the methods used in reducing these data from the time— 
histoiy records to the parameters presented in this paper is given in the 
appendix. 

The Reynolds number and the Mach number ranges of the models and 
the representative full-scale airplane are shown in figure 11. 


RESULTS AND DISCUSSION 


AERODYNAMIC COEFFICIENTS AND STAB ILI TY DERIVATIVES 

Lift 


The lift data are presented in the form of lift— curve slope Cj^ 

for various Mach numbers (fig. 12) as obtained from two models of the 
same configuration but having different center— of— gravity locations and 
different weights. The range of angle of attack in which data were 
considered for determining C]^ was ±15°. The lift coefficient varied 
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linearly with angle of attack in this range. As indicated in figure 12, 
Ci^ increased approximately 25 percent from the lowest test Mach num- 
ber (M = 0.88) to a Mach number of 1.00 and then decreased approximately 
15 percent from M = 1.00 to M = 1.20. The increase in lift-curve 
slope in going through the transonic region was evident for both of the 
models . Unpublished data obtained from wind-tunnel tests of a similar 
model for both high subsonic and low supersonic velocities have also been 
plotted in figure 12. 


Trim Lift Coefficient 


The variation of trim lift coefficient . with Mach number 

at different elevator deflections for model 1 is shown in figure 13(a) 
and for models 2 and 3, in figure 13(h). Different elevator settings 
for models 2 and 3 confirmed the as-sumption that Cj^. irn varied linearly 

with elevator deflection. These plots show an inherent characteristic 
of the model configuration to trim at negative lift coefficients between 
Mach numbers of 0.90 and 1.08. This was due to a basic untrimmed 
pitching-moment coefficient for the airplane at zero angle of 

attack and zero elevator deflection. Figure 14 shows a plot of . C^ as 

a function of Mach number. The asymmetry of the model configuration due 
to the vertical tail and the upswept rear of the body would indicate an 
expected positive Cm^ which was not in accord with test results. This 

negative trend of figure 14 however does agree with the data of 
reference 4. 


Change of Trim Lift Coefficient with Eespect 
to Elevator Deflection 

The rate of change in trim lift coefficient with respect to elevator 
deflection C]u is shown as a function of Mach number in figure 15 

°trim 

for models 1, 2, and 3. As would be expected, the values of Ct 

^trim 

for model 1 with the center of gravity at 25 percent mean aerodynamic 
chord were larger than those of models 2 and 3 with the center of gravity 
at 20 percent mean aerodynamic chord. Within the Mach number range 

covered by the tests, C]u remained fairly constant up to M = 0.86 

°trim 

at which point an abrupt reduction from 0.049 to 0.029 occurred 
between M = 0.86 and M = 1.00. A further decrease from 0.029 to 0.015 
occurred in Ci_ between M = 1.00 and M = 1.28. 

™trim 
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Hinge-Moment Coefficients 


The variations of hinge-moment coefficients with angle of attack 
and with elevator deflection are "both plotted as functions of Mach number 
in figure 1 6. The data points shown are for models 1 and 3* The theo- 
retical points shown on the plot -were calculated for a constant-chord 
partial— span control surface on a thin triangular wing as described in 
reference 5 • 

Calculations were made to determine the effect of elevon inertia on 
the hinge-moment coefficients. An extreme case showed the magnitude of 
the error to be negligible. Therefore, no such correction was applied 
to the data. Corrections were applied to eliminate the effect of phase 
lag between the hinge-moment coefficient and angle— of -attack curves and 
the effect of oscillations in elevon deflection due to angle— of -attack 
changes. Hinge-moment coefficients plotted as functions of angle of 
attack at a constant Mach number indicated that the variation was linear 
in the range covered by the tests (a = ±15°). 

Figure 1 6 shows that Cj^ increases from — 0.008 at M = O.85 

to -0.024 at M = 1.20. A corresponding increase from -0.015 to —0.037 
is shown for Ch§ between M = 0.85 and M = 1.05. Both curves indicate 

a gradual decrease in the low supersonic region. 


The value of the basic hinge-moment coefficient at zero angle of 

attack and elevator deflection Cy, is shown as a function of Mach 

0 


number in figure 17* 


The basic hinge-moment coefficient 


0u shows a 
~hQ 


reversal from positive to negative values at M = 0.95 and a tendency 
in the low supersonic region to return to positive values. The varia- 
tion of hinge-moment coefficient with elevator deflection was assumed 

to be linear in the solution of Cy, . 

n 0 


■ 1 Control Effectiveness 

A characteristic of the elevator used on the models can be seen in 
the plot of change in lift coefficient per degree of elevator deflec- 
tion Cj^ as a function of Mach number (fig. l8) . The parameter 

reaches a value of 0.022 at a Mach number of 0.96 and decreases to a 
value of 0.010 at M = 1.17, a reduction of about 55 percent through 
this speed range. Values of C^ are contained in' reference 4 and show 

good agreement with the flight— test values obtained in the high subsonic 
and low supersonic regions. 
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Two more parameters of longitudinal-control effectiveness for this 
configuration are shown in figures 19 and 20, change in trim angle of 


( X§r) , and change in 
AS /trim 


attack per degree of elevator deflection 
pitching-moment coefficient per degree of elevator deflection 

both shown as functions of Mach number. These two plots indicate an 
abrupt decrease in control effectiveness of the elevon between M = 0.90 
and M = 1.00. This reduction is of the order of 25 percent for 

^AS/trim* 

indicate a further gradual decrease in longitudinal-control effective- 
ness to M = 1.28, the highest Mach number reached by the flight tests. 
Values of were determined for the angle-of-attack r an ge 

between 10° and —8° . 


Cms 


and 35 percent for 


Above a Mach number of 1.00 the curves 


The effect of center-of-gravity location is apparent in both plots 
by the relative displacement of results obtained from model 1 with the 
center of gravity at 25 percent mean aerodynamic chord and models 2 
and 3 with the center of gravity at 20 percent mean aerodynamic chord. 
The more rearward location of the center of gravity reduced the value 

of Pm— and increased the magnitude of ( 

^ \ A6 /trim 


Longitudinal Stability 

Static stability .— When the controls are moved up and down in a 
square-wave type of motion, corresponding changes are produced in angle 
of attack and normal acceleration. The stability of the configuration 
is indicated by the period and the rate of decay of the short-period 
longitudinal oscillation when the controls are held fixed between pulses. 

The values of the period of the short-period oscillation induced 
by this abrupt control movement as determined from the time-history 
records are presented in figure 21 to show the variation of the period 
with Mach number for the models. The period decreased, indicating a 
stability increase, from a Mach number of 0.75, the lower test limit, 
to approximately M = 0.95. Above this speed the period continued to 
decrease but at. a much more gradual rate up to M = 1.28, the upper 
limit of the speed range covered by the flight tests. The period for 
model 1 was greater than that for models 2 and 3 throughout it's test 
range as would be expected since the center of gravity of model 1 
was 5 percent of the mean aerodynamic chord behind the center-of-gravity 
location for models 2 and 3. 
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The static— longitudinal-stability parameter in the form of the 
change in pitching-moment coefficient with respect to a change in angle 
of attack Cm is shown as a function of Mach number in figure 22 


a 


for values between ±0.30. The determination of involved the 


use of the period of the short— period oscillations as a primary factor. 
The value of increased from a minimum of -0.0095 at M = 0.85 to a' 

rnavimum of -0.0l62 at M = 1.15 for models 2 and 3 with the center of 
gravity at 20 percent mean aerodynamic chord. An investigation of the 
change in Cn^ due to a 5 -i percent change in center-of— gravity location 

shows that C m for model 1 is lower than would be expected from a com- 


parison with models 2 and 3* Data concerning the evaluation of 

were carefully rechecked and indications are that the seemingly low 
values of Cm^ were due to accumulative errors within the accuracy of 

determining the physical characteristics used to calculate this parameter. 


Figure 23, a plot of aerodynamic-center position against Mach 
number, also indicates the variation of the static longitudinal stability 
of this configuration. The aerodynamic-center positions for model 1, 
however, were 2^- percent of the mean aerodynamic chord ahead of 
models 2 and 3 . The more forward aerodynamic-center locations for 
model 1 were a result of the low values of obtained for this model. 

This difference, however, is within the accuracy of aerodynamic-center 
location usually obtained from flight and wind-tunnel data. 


The three parameters discussed in the preceding paragraphs (period, 
0m , and aerodynamic— center position) show that the static longitudinal 


stability of this configuration increased through the transonic region 
from a minimum at about M = 0.82 to a maximum value at M = 1.15* 


Dynamic stability .— A qualitative evaluation of the dynamic stability 
may be made by inspection of the damping of the short— period oscillation 
induced by the abrupt control movement. Damping is represented by the 
parameter T \j2> the time required to damp to one-half amplitude, and is 

presented in figure 2k as it varies with Mach number. Since the flight- 
test models were not dynamic— scale models, the results presented in 
figure 24 are applicable to the full— scale airplane only after corrections 
are applied. Models 2 and 3 with the center of gravity at 20 percent 
mean aerodynamic chord showed more rapid damping characteristics than 
model 1 with its center of gravity at 25 percent mean aerodynamic chord. 
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The total damping factor Cm- - + Cm. _, which is a measure of the 

9 c ooc 

2Y 2V 

dynamic stability of the configuration expressed nondimen sionally, is 
presented in figure 25 as a function of Mach number. Model 1 with. the 
more rearward center— of— gravity location indicated less tendency to 
damp throughout the flight— test speed range than did models 2 and 3 • 

It will be noted that there is considerable scatter in the damping 
data. This type of scatter may also be expected for full-scale airplane 
conditions inasmuch as the present data were obtained in free flight and 
a.l 1 the aerodynamic factors that affect damping were properly integrated 
into the motion of the models. 


Directional Stability 

Only models 1 and 2 were instrumented to obtain transverse acceler- 
ations. Model 2 apparently had some directional asymmetry causing it 
to develop a small positive side force throughout the flight. This 
effect approximately doubled at Mach numbers below 0.90. Model 1 did 
not exhibit any such consistent side— force variation, the side forces 
on model 1 resulting from an occasional disturbance. Neither model 
showed divergence nor continuous oscillations thus indicating positive 
directional stability. 


FIXING QUALITIES 

Longitudinal Trim Characteristics 


Trim angle of attack .— The angle- of attack for trimmed level flight 
required for this configuration is presented as a function of Mach number 
in figure 2 6. Curves are shown for center— of— gravity locations at 20 
and 25 percent of the mean aerodynamic chord for both sea— level flight 
and flight at an altitude of 40,000 feet. 

The trim angle of attack shows a consistent small decrease with 
increasing speed except in the region between M = 0.90 and M = 0.95* 

Control position for trim .— The characteristics of the elevator 
control in level flight are presented in figure 27 in the form of the 
variation of the elevator position required for trim with Mach number. 
Control— position trim change is manifested between a Mach number of 0.87 
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and 0.95 at sea level and 40,000 feet. The control— position trim 
change is a function of variation of out— of— trim pitching moment with 
Mach number, change in control effectiveness, and movement of the 
neutral point. The resultant change in trim, a tucking-ainder tendency, 
appears to be of moderate magnitude. For example, at 40,000 feet a 
maximum up-elevator angle of about 5° is required for trim at a Mach 
number of 0.95* 

An evaluation of the stick— fixed maneuver point in the Mach number 
range between 0.80 and 1.20 indicated that the point is well behind the 
most, rearward center— of— gravity position and the requirements are met 
for maneuvering stability in reference 2. 

Figures 28, 29, 30, and 31 present the control positions required 
for maneuvering at various accelerations as functions of Mach number. 
Figures 28 and 29 are curves for sea level and 40,000 feet, respectively 
with the center of gravity located at 25 percent mean aerodynamic chord. 
Figures 30 and 31 are also curves for sea level and 40,000 feet, respec- 
tively, but with the center of gravity located at 20 percent mean aero- 
dynamic chord. For example, figure 29 shows that at a Mach number 
of 1.20 an up-elevator deflection of 10° would produce 5g normal 
acceleration at 40,000 feet with the center of gravity located at 
25 percent mean aerodynamic chord. 

Longitudinal control forces .— The elevator control force required 
for trim in straight and level flight at various Mach numbers is pre- 
sented in figure 32 for sea— level flight and in figure 33 for flight 
at 40,000 feet. The stick force per g is presented in figure 34 as a 
function of Mach number. The stick forces are based on a conventional 
airplane configuration with 2° of elevator deflection for 1 inch of 
stick movement. Therefore, these data indicate the power required of a 
control-boost system with no balancing and trimming devices. For 
example, with the center of gravity at 25 percent mean aerodynamic chord 
at a Mach number of 1.20 the stick force per g based on measured hinge 
moments is about 900 pounds per g. 

The variation of elevator control force for trim with Mach number 
(fig. 32) indicated that pull forces were required at all speeds below 
the trim speed and push forces required at all speeds above the trim 
speed within the range of Mach numbers from 0.95 to 1.20. The opposite 
is true for Mach numbers from 0.80 to 0.95, but the elevator angle for 
trim in this range, of Mach number increases with increasing Mach number. 
Thus the stick force would be in the correct sense with respect to stick 
movement throughout the transonic region. 

The elevator hinge-moment data obtained for model 1 indicate a 
force reversal at high angles of attack (a s? 15°) at Mach numbers 
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"below 0.90. Model 2, which flew at angles of attack of about 7° 
at M = 0.90, did not show a hinge-moment reversal but did indicate 
hinge moments near zero. 


Longitudinal Control Effectiveness 

The variation with Mach number of the' normal acceleration produced 

per unit of . elevator deflection (^) is presented in figure 35* At 

NAb/trim 

sea level a large variation in elevator effectiveness was apparent from 
subsonic to low supersonic speeds with minimum effectiveness occurring 
at a Mach number of 1.06 for model 1 with the center— of— gravity location 
at 25 percent mean aerodynamic chord and at a Mach number of O .98 for 
models 2 and 3 with the centei>-of— gravity location at 20 percent mean 
aerodynamic chord. Sufficient control for maneuvering is available as 
indicated by the fact that 10° elevator deflection will produce 5g acceler- 
ation at a Mach number of 1.20 at 40,000 feet with the center of gravity 
located at 25 ^ percent mean aerodynamic chord. 


Dynamic Stability 

The characteristics of the stick— fixed short— period longitudinal 
oscillations are presented in figures 36 to 38 for a full-scale con- 
figuration. U. S. Air Force specifications for stability— end— control 
characteristics of airplanes (reference 2) require that the short— period 
dynamic oscillation of normal acceleration produced by moving and quickly 
releasing the elevator shall be damped to l/lO amplitude in one cycle 
(based on free controls). The damping characteristics for the full-scale 
configuration have been evaluated for the control— fixed condition although 
there is a slight oscillation in the control position due to hinge-moment 
effect which is apparent in figures 8 to 10. The fixed-control charac- 
teristics would probably dictate the behavior of this airplane since 
it would require some kind of control— boost system to aid the pilot in 
overcoming the extremely large stick forces encountered in maneuvering. 
Figure 36 , which gives the cycles required to damp to l/lO amplitude 
as a function of Mach number at sea level and 40,000-foot altitude, 
indicates that this tailless design -would more than meet such a require- 
ment at sea level for the speed range tested since the longitudinal 
short— period oscillation will damp to an average value of l/l6 amplitude 
in one cycle. At 40, 000. feet, however, the damping does not meet the 
requirement at any speed in the range tested. 

Figure 37 shows the time required for the short— period oscillations 
of the full-scale airplane to damp to l/2 amplitude as a function of 
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Mach number at sea level and 40,000 feet. As can be seen from this 
figure, Tjy 2 decreases through the transonic region and reaches a 

relatively constant value at about M = 1.20. Period as a function of 
Mach number is shown in figure 38 for sea level and 40,000 feet. 

Both T i/ 2 a and- period indicate increasing stability for the configu- 
ration with increasing Mach number in the transonic and low supersonic 
speed range . 


CONCLUSIONS 


From the results of a flight investigation made to evaluate the 
aerodynamic characteristics and flying qualities of models of a tailless 
triangular-wing airplane configuration, the following general conclusions 
are indicated for the Mach number range between O .75 and 1.28. 


AERODYNAMIC CHARACTERISTICS 


1. The lift coefficients varied, linearly in the angle— of -attack 
test range of ± 15 °. The lift-curve slope Cj^ varied from 6.0 at 

a Mach number of 0.88 to a maximum of 0.055 at a Mach number M of 1.0 
and then decreased to O.OV 75 at a Mach number of 1.20. 

2. The hinge-moment coefficient per degree of angle of attack 
increased 200 percent between M = O .85 and M = 1.20; whereas the 
hinge-moment coefficient per degree of elevator showed a corresponding 
rise of 150 percent between M = 0.85 and M = 1.05. Both of these 
values showed a gradual decrease in the -low supersonic region. 

3. The elevator effectiveness decreased by approximately ho percent 
from. a Mach number of 0.9 to 1.25 . For example, with the center of 
gravity at 20 percent mean aerodynamic chord, the rate of change of 
pitching-moment coefficient with elevator deflection Cj^ at a Mach 

number of 0.9 was — 0.015 and at a Mach number of 1.25 was -0.009. 

4. The configuration tested possessed static longitudinal stability 
throughout the Mach number range covered by these flight tests. The 
value of Cjjj^ (rate of change of pitching-moment coefficient with angle 

of attack) increased from a minimum at M = 0.80 to a max imum 

^ — 1*15 with the center of gravity at 20 percent mean aerodynamic 
chord. 
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5- The aerodynamic center moved very gradually from a minlTmiTn 
of 42 percent of the mean aerodynamic chord at a Mach n umb er of 0.80 to 
a max imu m of 5 4 percent of- the mean aerodynamic chord at a Mach number 
of 1.15. 

6. The damping parameters and coefficients indicated that the con- 
figuration possessed dynamic longitudinal stability throughout the test 
speed range. 

7. The models exhibited directional stability throughout the angle- 
of— attack and speed ranges of the tests . 


FLYING QUALITIES 


1. There is ample control for trim in level flight at sea level 

or at altitude. At 40,000 feet a maximum up-elevator angle of about 5° 
is required for trim at a Mach number of 0.96. The transonic trim 
change, a tucking-amder tendency, appears to be mild. 

2. The elevator control remains effective in changing lift or 
angle of attack over the entire speed range. The effectiveness of the 
elevator in changing angle, of attack, however, is reduced to about half 
of its subsonic value at supersonic speeds. This change of effective- 
ness occurs gradually. 

3. With the center of gravity at 25 percent mean aerodynamic chord 
the normal acceleration produced per degree elevator deflection is such 
that about 10° up-elevator deflection is required to produce a 5g 
acceleration at 40,000 feet at a Mach number of 1.2. The corresponding 
stick force per g based on the measured hinge moments is about 900 po und s 
per g, a figure which gives an indication of the power required of a 
control— boost system. 
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4. The damping of the short— period longitudinal oscillation is 
adequate over the speed range for the sea— level condition (of the order 
of 1 cycle to damp- to l/l6 amplitude). At b0,000 feet, however, the 
damping is inadequate in respect to the U. S'. Air Force requirements. 


Langley Aeronautical Laboratory . 

National Advisory Committee for Aeronautics 
Langley Air Force Base, Va. 
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APPENDIX 

REDUCTION OF DATA AND METHODS 
OFANALTSIS 
Reduction of Data 


Mach number .— The total pressures obtained from the telemeter 
records were reduced to Mach number by use of the following equations: 

Subsonic 


H 

P 




M‘ 


i 


2L 

7-1 


( 1 ) 


Supersonic 



( 2 ) 


where p, free— stream static pressure, was obtained from the reference 
static— pressure record in conjunction with radiosonde data. Models 1 
and 3 reached a maximum altitude of 4000 feet while model 2 attained a 
maximum of 4-700 feet. The Doppler velocimeter radar unit served as an 
independent check of the Mach number obtained by reduction of the total 
and reference static pressures. 

Angle of attack .- Since angle-of -attack data were measured at a 
point some distance ahead of the center— of— gravity location, it was 
necessary to correct these data for flight— path curvature and angular 
velocity as described in reference 1. The following equation was used: 


a = 


+ | ( 1844 — — + 

v \ s v 


dt / 


(3) 
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Control position .— Prior to the flight test of each model a static 
hinge-moment calibration of the control system was conducted to determine 
the amount of twist that would be encountered in the elevons and control 
linkage under aerodynamic loads. The elevons were loaded at two spanwlse 
stations and readings were taken at five points to measure the amount 
of twist or deflection induced. Control— position data recorded during 
flight were corrected by the calibration obtained from the static test. 


Analysis of Aerodynamic Coefficients 
and Derivatives 

Basis .— The methods of analysis used herein apply to the free 
oscillation resulting from a step function disturbance. The disturbance 
was created by an approximately square-^wave type of moti.on of the eleva- 
tors moved abruptly between limit stops. The complete derivation of the 
equations used will not be given herein but the basic equations of motion 
are as follows: 


Vm(0 - a) = + %5^57.3qS 

I Y 0 = + C^a + C^0 + 0^57.3(135 


W 

(5) 


In order to simplify the analysis and to permit the determination 
of equations for the more important aerodynamic derivatives, a number of 
assumptions are necessary. It is assumed that during the time interval 
over which each calculation is made the following conditions exist: 

The forward velocity is constant and the aerodynamic forces and moments 
vary linearly with the variables a, 8, and 0. 

Figure 39 is a schematic plot showing a typical record of the con- 
trol position and lift-coefficient responses- following step deflections 
of the aircraft control surfaces. At least three complete peaks of each 
disturbance were necessary to obtain the trim lines shown in the oscil- 
lations. Where three complete peaks were not present, sufficient 
accuracy in placing the trim line could not be ascertained and such data 
carried little weight in the analysis. 

Lift-curve slope .— Several methods were tried for determining the 
lift-curve slope with respect to angle of attack. The most expeditious 
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method found was to measure the instantaneous slopes 


— — and ^ at 
dt dt 


a given Mach number. Care was exercised in using only the portions of 
the lift coefficient and angle— of-attack time-history curves where 
slopes could be accurately ascertained. The effect of lift due to the 


flexibility of the elevator was eliminated by correcting 


d^ 

dt 


for the 


lift due to the deviation of the elevator deflection from a fixed value 
at an instantaneous time. The following relation exists: 



AC I 2 _ 1 “ C I 6 a5 2-1 
** 2-1 


(6) 


where AC^ is the change in between Cp^ and taken over 

a relatively straight portion of the lift time history as indicated in 
figure 39 and A&2— i and ^2—1 8X6 incremental changes in 5 and a 

over the same time interval. The value of Cp^. used in this equation 

was a first approximation. Successive approximations and evaluations 
were unnecessary. 


After determining the corrected value of it was then possible 


a 


to determine an exact value for C j^ , the lift— curve slope due to the 


elevons , from thq portions of the time histories where the controls were 
moving from one extreme position to the other. The following relation 
exists: 


(^-lXrim CL a( a2 O trim 

^ = (Vl) trim 


* (7) 


The values of a, 5, and C L are trim values as illustrated in figure 39. 

Although / oc _ \ is not included in the illustration, it would be 

\ 2 ■ l/trim 

obtained in the same manner. 
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The variation of trim lift coefficient with respect to elevator 


deflection C T 


was found by the same method used to find C- 


■’trim 


V 


The equation is 


’trim 


K-i) 


trim 


(8) 


trim 


and the quantities used are again illustrated in figure 39. 

The trim lift coefficients C T corresponding to the trim 

-^rim 

elevator deflections encountered in the tests were plotted against Mach 
number in figure 13* Trim lift coefficients for elevator settings 
between ±5° were derived by using a linear relation between lift coef- 
ficient and elevon deflection at a constant Mach number. 


The basic untrimmed pitching-moment coefficient was calculated 

from the conventional moment-coefficient equation solved for 
follows: 


as 


C ™0 “ “Sx^rim (^^g^trim (9) 


Cm^ and C mr- were considered linear in the range of the tests. The 
second term was eliminated by taking values of . for zero degrees 

of elevator deflection and dividing the first term by Cj, to mak e 
a function of the trim lift coefficient, or 


and 


Cffl 0 _Cm a . a trim &=0 

S = “ C^( CL trim) 5=0 


(10) 


( 11 ) 


The values of 0m^ were obtained as described in a following section. 
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Hinge momenta .— Hinge-moment data were reduced to coefficient form 
and plotted directly against angle of attack for both up and down elevator 
to obtain an approximate value of . This value was used in conjunc- 


tion with the change in 5 . due to changes in a to correct the values 
of the' total hinge moment for constant elevator deflection as follows: 


( H= e = 0,1 + 


( 12 ) 


The values obtained were plotted as functions of time as were the values 
of angle of attack. A method was derived to eliminate the effect of 
phase lag between the two variables. Constant values of were chosen 

on each side of an oscillation peak and a mean value of a corresponding 
to the constant value of C^ was determined analytically and graphically. 

Finally the corrected values of C^ and a were plotted for up and down 
elevator to determine C h ^ and C h& . Indications were that these values 

were linear and C^q, the hinge-moment coefficient at zero angle of attack 
and elevator , was determined by direct interpolation. 

Control effectiveness .— The variation of trim angle of attack with 

elevator deflection ( && \ was found by using the method illustrated 

VA5/ trim 

by figure 39 and the following equation 


'^8 /trim 




(13) 


The resulting values were used to obtain the control-effectiveness 
parameter Cm^ at a constant angle of attack in the manner given as 

follows: 


KLk = 


(14) 


The solution of Cm^ is presented in the discussion of longitudi nal 
stability. which follows. 
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Longitudinal stability .- Evaluations of the longitud inal stability 
were obtained by analysis of the short-period oscillations induced by 
the abrupt control movements and shown in the angle-of -attack curves in 

the time histories. The solution of , the static longitudinal sta— 

a 

bility. derivative, is obtained from the following equation as derived 
from the simultaneous solution of the two equations of motion: 



flexibility and the second-order effects from the two-degree s-of- freedom 
method of analysis were neglected since they constituted less than 
0.5 percent of the results. 


The periods of the short-period oscillation P were read from the 
time-history curves and the time to damp the amplitudes to one-half 
magnitude was determined by the use of the following formula: 


, = 0.693 P 

1/2 2 log e (a-j/A^ 


(16) 


where A^ and Ag were successive amplitudes above and below the 
neutral axis of the angle-of -attack time history at the point where T-jy 
was sought. 

The quantities CL and C j^ , corrected for the effect of elevator 

oscillations, were used in conjunction with the model center-of-gravity 
locations to determine the aerodynamic-center positions in percentages 
of the mean aerodynamic chord. Eeference 6 presents the following 
relation: 


Aerodynamic center = Center of 


gravity - 



(17) 
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The dynarni c— longitudinal— stab ility data were reduced to the form 

of Cni._ + C m by the following equation derived from the simultaneous 
6c ac 

2V 2V 

solution of the two equations of motion: 


2V 



8l I 0.693 

P T 

pVSc I 1/2 


57.3Cl^pVB 

4m 


(18) 


Flying-Qualities •■Analysis 

Variation with Mach number of the control position required for 

trim in level flight.— The trim lift coefficient Cl for 0° elevator 

: trim 

deflection was obtained by plotting values of . corresponding to con- 
stant positive and negative elevator deflections against Mach number, and 
the variations were considered to be linear in the test range. These 
values were taken from the trim values of and 6 obtained from the 

time-history data of the flight tests of the three models. The value 

of Cl for 6=0° was obtained by interpolation. Values of Cl 

trim 

for level flight for the full-scale airplane were obtained from the 
relation (Cl)^ = — The difference between (Cl)^ for straight and 

level flight and Ct for 8 = 0° was divided by Ct to give 8 

trim °trim ■ 

for straight and level flight for various Mach numbers. 


8 = 


( C;L )lg ( CL trim) d= 


8 = 0 ° 




trim 


(19) 
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Elevator control force for trim against Mach number .— A value of 
deflection of elevator per inch of stick movement for a high-speed 
fightei>-type airplane was assumed to "be 


— =2° per inch 
x 


Values of hinge moment were obtained from the time-history plots of models 
for corresponding values against Mach number. The method for 

determining trim lines was the same as in figure 39. The value of ( — ) 

VAB/ trim 

was obtained from 


(iS\ = ( Hg Bl )trlm 

VAS /trim (b 0 — 5. \ 

V 2 x / trim 


At a given Mach number a value of hinge moment was read at a given elevator 


deflection and corrected to the 6. 


trim 


for straight and level flight at 


sea— level conditions by 


\rin = Ml . ' 

If the hinge moment for S trlm for straight and level flight at sea- 
level conditions is known, the elevator control force is obtained by 



H 5 
57-3 x 


( 22 ) 


where H has heen corrected to full scale 
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Change in normal acceleration for a corresponding change in elevator 

deflection ( — ] against Mach number.— The values of Ct for level 

\^/trlm 

flight for various Mach numbers were divided by Ct so that for lg 

^trim 


££> = 


H, 

0l ®trlm 


(23) 


Stick force per g against Mach number .— The change in elevator 
deflection required for a change in normal acceleration of lg, reciprocal 

of ( — ) , was multiplied by — — — to obtain the change in hinge 

AB ' tr im ^tr im 

moment required for a change in normal acceleration of lg. Then 
AF 

for — in pounds per g 

g 


AF = /A)\ AH 6 1 

8 V 8 /trim ABtrim x 57-3 


(24) 


Dynamic stability .— The dynamic stability of the airplane in terms 
of period and damping of the short— period longitudinal oscillations was 
determined from the oscillations of the model corrected to full-scale 
conditions. 

The correction factors were determined from a two— degree— of— freedom 
method of analysis of the motion which assumes no changes in forward speed 
during the oscillation. The period of the oscillation for the airplane 
in terms of period for the model was obtained from a ratio of the C m^ 
equations for the two as 



'^a ^m c m Pm 
Pa 


S a c a % 


m 


( 25 ) 
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The time to damp to one-half amplitude for the airplane was determined 
by the following relationship: 


Cm.- + Qm.. — = 
Sc ac 


-BTy 


2Y 


^ paMSc 2 


0-693 

r i/2 


57.3C Itt PaMS 

4m 


(26) 


and equated for model and airplane as follows: 


0.693 

T l/2. 


57.3C La P a a a MS a 


- 2 
c c 




L Y^a 
m 


“a 


% % c m 2 

a 


% m a a p a S a c a 

! % ^Pm^m^m^ 
a 


fM>3\ ( 27 ) 


Flying-qualities specifications require that the short-period 
oscillations damp to l/lO amplitude in one complete cycle. This value was. 
determined from the relation 


C l/lO a 


3.325T. 


1/2* 


(28) 


for the representative full-scale airplane. 
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TABLE II 

WEIGHT AND BALANCE DATA FOE MODELS 1, 2, AND 3 AND FULL-SCALE 
TEIANGULAE-WING, TAILLESS CONEIGUEATION 


Model 


Weight 

(lb) 



Wing loading 
(lh/sq ft) 

Center— of -gravity 
position 
(percent M.A.C. ) 

Moment of inertia 
(slug-ft 2 ) 

X Y 

Socket fuel 

included in the models 

30.1 

29.7 

17-52 

30.4 

24.0) 

17.89 

29.9 

24.1) 

18.76 * 


Models without rocket fuel 


1 

2 

3 

182.50 

184.25 

181.25 

29.2 

29.5 

29.0 

25.0 ' 

20.0 

20.0 

16.65 

17.10 

17.93 

Normal gross 
weight (lb) 

Wing loading 
(lb/sq ft) 

Center-of-gravity 
position 
(percent M.A.C.) 

Moment of inertia 1 

( 8lug— f t 2 ) 


Full-scale configuration 


11,600 

11,600 


27.3 

25.0 

27.3 

20.0 






















-Elevon hinge 
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Figure 2.- Side view of model 
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Figure 3»- Bottom view of model. 
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L- 57682 


Figure 4 .- Model-booster combination on launcher. 
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PICK-UP OSCILLATORS 


NORMAL 

ACCELEROMETER PICKUP 


LONGITUDINAL 
ACCELEROMETER PICKUP 




Sf 

AIR PRESSURE 
SWITCHING MOTOR 

AIR SERVO 


ELEVON HINGE-MOMENT 
PICKUP 


ELEVON ACTUATING ARM 
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TRANSMITTER — 

MODULATOR . 
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TOTAL HEAD PRESSURE 
PICKUP 
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BOX 


TRANSVERSE 
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CONTROL POSITION 
PICKUP 


Figure 6.- Instrumentation arrangement on a typical model 




Figure 7*- Side view of model showing bent angle-of-attack sting. 
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Mach number, M 


Figure 11,- Variation of Reynolds number with Mach number for flight of 
the full-scale airplane at two altitudes and for flight model test 
data. 
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Figure 12.- Change in lift coefficient with respect to angle of attack as a function of Mach number. 
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Mach number, M 


Figure lU.- Variation of the basic untrimmed pitching-moment coefficient 

with Mach number. 
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Mach number, M 


Figure 15*- Change in trim lift coefficient with respect to elevator 
deflection as a function of Mach number. 
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Figure l 6 .~ Variation with Mach number of the change in hinge-moment 
coefficients with respect to angle of attack and elevator 
deflection. 



NACA RM L 9 LO 7 


53 



Mach number, M 


Figure 17.- Variation of the basic hinge-moment coefficient at zero 
angle of attack and zero elevator deflection with Mach number. 



Mach number, M 


Figure 18.- Cheuige in lift coefficient with respect to elevator 
deflection as a function of Mach number. 
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Figure 19*- Change in trim angle of attack with respect to elevator 
deflection as a function of Mach number. 



Mach number, M 


Figure 20.- Change in pitching -moment coefficient at a constant angle 
of attack with respect to elevator deflection as a function of 
Mach number. 
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Figure 21.- Period of the short-period longitudinal oscillations as a 

function of Mach number. 
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figure 22.- Pitching-moment-coefficient slope with respect to angle of 
attack as a function of Mach number. 
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Figure 23.- Aerodynamic-center position in percent of the mean 
aerodynamic chord as a function of Mach number. 
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Figure 2k.- Variation with Mach number of the time required for the 
short-period longitudinal oscillations to damp to one-half 
amplitude . 
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Figure 27.- Variation with Mach number of the elevator deflection 
required for trim in level flight at different altitudes and 
center-of-gravity locations . 



Mach number, M 


Figure 28.- Variation with Mach number of the elevator deflection 
required to produce various maneuvers at sea level with the 
center of gravity at 25 percent mean aerodynamic chord. 
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Figure 29-- Variation with Mach number of the elevator deflection 
required to produce various maneuvers at 40,000 feet with the 
center of gravity at 25 percent mean aerodynamic chord. 
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Figure 31- - Variation with Mach number of the elevator deflection 
required to produce various maneuvers at 1+0,000 feet with the 
center of gravity at 20 percent mean aerodynamic chord. 
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Figure 32.- Variation with Mach number of the elevator-control force 
required for trim in level flight at sea level for two center-of- 
gravity locations. 



Figure 33 . - Variation with Mach' number of the elevator-control force 
required for trim in level flight at 1+0,000 feet for two center- 
of-gravity locations. 
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Mach number, M 


Figure 3^>- Variation of the stick force per g with Mach number for 
different center-of -gravity locations. 



Figure 35*- Variation with Mach number of the normal acceleration 
produced per -unit elevator deflection at different altitudes and 
center-of-gravity locations. 
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Mach number, M 


Figure 36.- Variation with Mach number of the cycles required for the 
short-period longitudinal oscillations of the full-scale configu- 
ration to damp to one -tenth amplitude at different altitudes and 
center-of-gravity locations. 
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Figure 37*- Variation with Mach number of the time required for the 
short-period longitudinal oscillations of the full-scale configu 
ration to damp to one -half amplitude at different altitudes and 
center-of-gravity locations. 
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Mach number, M 


Figure 38.- Variation with Mach number of the period of the short- 
period longitudinal oscillations for the full-scale configuration 
at different altitudes and center-of-gravity locations. 




Figure 39 -- Typical section of time history as used for data analysis 






